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Abstract. The kinetic decoupling of weakly interacting massive particles (WIMPs) 
in the early universe sets a scale that can directly be translated into a small-scale 
cutoff in the spectrum of matter density fluctuations. The formalism presented here 
allows a precise description of the decoupling process and thus the determination of this 
scale to a high accuracy from the details of the underlying WIMP microphysics. With 
decoupling temperatures of several MeV to a few GeV, the smallest protohalos to be 
formed range between 10~^^ and almost 10~^ solar masses - a somewhat smaller range 
than what was found earlier using order-of-magnitude estimates for the decoupling 
temperature; for a given WIMP model, the actual cutoff mass is typically about a 
factor of 10 greater than derived in that way, though in some cases the difference may 
be as large as a factor of several 100. Observational consequences and prospects to 
probe this small-scale cutoff, which would provide a fascinating new window into the 
particle nature of dark matter, arc discussed. 
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1. Introduction 

The combined data of cosmic microwave observations [1], distance measurements of type 
la supernovae [2] and the baryon acoustic oscillations in the distribution of galaxies 
[3], along with a wealth of other observations, provide ample evidence for a minimal, 
6-parameter cosmological ACDM model, where the dominant contribution to the total 
energy density of the universe today are dark energy and cold dark matter (CDM). While 
the nature of the latter (and even more so of the former) remains unknown, weakly 
interacting massive particles (WIMPs) provide a class of particularly well motivated 
CDM candidates that, thermally produced in the early universe, naturally acquire the 
necessary relic density today (for reviews on particle DM, see [4, 5]). It is interesting to 
note that WIMPs automatically arise as a by-product of almost any attempt to address 
the problems of the standard model of particle physics (SM) that become apparent at 
the TeV scale, viz. its extreme UV sensitivity in the scalar sector, a scale that soon 
will be accessible for the first time by the LHC The most popular, and arguably most 
well- motivated scenario in this context is supersymmetry [4], but interesting alternatives 
include, e.g., scenarios with large extra dimensions [6] or little-Higgs models [7]. For 
concreteness, the focus will often be on the supersymmetric neutralino, but the main 
results - such as the mass range of the smallest WIMP protohalos - will be equally 
applicable to, and include, other WIMP DM candidates. 

After the WIMP number density has frozen out in the early universe, setting the 
DM relic density today, scattering events with the much more abundant SM particles 
still keep the WIMPs extremely close to thermal equilibrium. The complete decoupling 
from the thermal bath happens only considerably later [8, 9], when the temperature has 
dropped by another factor of ~ 10 — 1000, from which moment on the DM particles 
can stream freely and thereby wash out any matter density contrasts on small scales 
[10, 11, 12, 13]. Acoustic oscillations during and after kinetic decoupling also give rise 
to a damping of the power spectrum of matter density fluctuations [14, 15]; whether or 
not this damping is stronger than the one induced by free streaming, as we will see later, 
depends on the details of the WIMP model. In any case is it the kinetic decoupling 
scale that determines the cutoff in the density spectrum and thus provides a direct 
link between the size of the smallest gravitationally bound structures and the CDM 
particle nature. To make proper use of this link, it is necessary to develop a detailed 
understanding of the decoupling process so as to allow a precise determination of the 
decoupling temperature. 

This article is organized as follows. Section 2 briefly introduces the supersymmetric 
neutralino as DM candidate, which for most of the article will be used as a prototype 
for the more general class of WIMP DM. In Section 3, a detailed description of the 
decoupling of WIMPs from the thermal bath is given, together with a prescription for 
a highly accurate determination of the kinetic decoupling temperature. The subsequent 
evolution of matter density fluctuations until the formation of the first protohalos is then 
discussed in Section 4. Observational consequences and prospects to probe the power 
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spectrum at the smallest scales are outlined in Section 5. Here, also the situation for 
other (WIMP and non-WIMP) DM candidates are briefly discussed. Section 6, finally, 
concludes. The two appendices provide analytic expressions for the collision term in the 
Boltzmann equation, describing the scattering between arbitrary DM candidates and 
SM heat bath particles, as well as the full scattering amplitude for neutralino-fermion 
scattering in supersymmetry. 

2. Supersymmetry and neutralino dark matter 

Major motivations to consider supersymmetric (SUSY) extensions to the SM include 
the (at least in the unbroken theory) perfect symmetry between fermions and bosons, 
which not only successfully addresses the unification of gauge couplings at the scale 
of grand unified theories (GUT) but also, for realistic schemes of SUSY breaking, the 
aforementioned hierarchy problem of the SM. The conservation of i?-parity, introduced 
to ensure the stability of the proton against decay, means that the lightest SUSY particle 
(LSP) has to be stable. In most models, the LSP is the lightest neutralino, a linear 
combination of the superpartners of the gauge and Higgs fields, 

X = Xi = NuB + NuW^ + NisH'i + N^.H^ , (1) 

which provides a perfect WIMP DM candidate (for a review, see [4]). 

A popular way of assessing the range of possible phenomenologies for neutralino 
DM is to study a 7 parameter version (at the electroweak scale) of the minimal 
supersymmetric extension of the SM (MSSM), where /i denotes the Higgsino mass 
parameter, tan/3 the ratio of vacuum expectation values of the two Higgs doublets. Mi, 
M2 and M3 the gaugino mass parameters, and rriA the neutral pseudoscalar Higgs mass. 
Further parameters are tuq, At and A^, which are defined through a simplifying ansatz 
to describe the soft trilinear couplings A and the soft sfermion masses M (which in 
general are 3x3 matrices in generation space): Mq = M[/ = Nl^ = '!^e = = "^o^, 
Au = diag(0, 0, At), Ad = diag(0, 0, A;,), A^; = 0. No CP-violating phases other than 
the CKM phase of the SM are allowed for; as a natural further simplification, the GUT 
condition for the gauge couplings, leading to Mi = | tan^ 6'^M2 ~ |M2 is used. 

A more restricted, but in some sense more natural, set of parameters defines 
models of minimal supergravity (mSUGRA), where gravity is supposed to mediate SUSY 
breaking. In these models, parameters are given at the GUT scale; the mass spectrum 
and couplings at the electroweak scale are then calculated by solving the renormalization 
group equations. The parameters of mSUGRA models are the universal scalar and 
gaugino masses, mo and mi/2, a common trilinear term Ao, tan/3, and the sign of /i. 
From a cosmological point of view, the regions in parameter space that are particularly 
interesting - as they correspond to models giving the correct relic density - are: The 
bulk region at low tuq and mi/2; the funnel region ~ 2m^, where annihilations 
in the early universe are enhanced by the presence of the near-resonant pseudoscalar 
Higgs boson; the hyperbolic branch or focus point region where mo ^ mi/2; the stau 



Particle DM and small-scale structure 



4 



coannihilation region where ^ ; and finally the stop coannihilation region (arising 
when Aq ^ 0) where m^^ ^ m^. 

To assess the range of decoupling scales for neutralino DM, an extensive scan over 
the MSSM and mSUGRA parameter space was performed (the same as described in 
[16]), resulting in about 4 ■ 10^ models that satisfy all current bounds from accelerator 
physics and at the same time feature a relic density within the 3(T bound of the combined 
WMAP, supernovae and baryonic acoustic oscillation data [1]: 

0.103 < f^^/i^ < 0.123. (2) 

For the numerical calculations in this article, the DarkSUSY code was used (see [17] for 
sign conventions and other details) which implements FeynHiggsFast [18] for the Higgs 
boson masses and decay widths, and Isajet [19] for the solution of the renormalization 
group equations (RGE) and for the mass spectra in mSUGRA models. The routines for 
the calculation of the decoupling scale as described in the next Section will be available 
with the next release of DarkSUSY [20]. 

3. Thermal decoupling of WIMPs 

At high temperatures in the early universe, T > m^, chemical equilibrium is maintained 
by the detailed balance between WIMPs annihilating into SM particles and the reverse 
process of SM particles creating WIMPs. In addition, local thermal equilibrium is 
guaranteed by highly frequent elastic scattering processes. Quantitatively, the evolution 
of the WIMP phase-space density / is described by the Boltzmann equation; denoting 
with = {E, p) the comoving WIMP momenta and assuming a fiat Friedmann- 
Robert son- Walker spacetime, it reads: 

E{d,-Hp-V^)f = C[f], (3) 

where H = a/a is the Hubble parameter and a the scale factor. The right-hand side is 
a source term, also known as the collision term, that describes any changes in / that 
cannot simply be attributed to the geometry of the expanding spacetime. 

Dividing Eq.(3) by E and integrating it over p, it becomes the familiar evolution 
equation for the WIMP number density n: 

dtn + ?,Hn = a'^dt {a^n) = -{av) {n^ - nl^) , (4) 

where Ueq denotes the number density in chemical equilibrium. Note that only non- 
number conserving contributions to C[f] survive the integration over p, which is the 
reason for the appearance of the thermally averaged annihilation rate (av) on the right- 
hand side. Incidentally, this equation takes the same form, after replacing (erf) with 
an effective annihilation rate, even when coannihilations with particles close in mass to 
are taken into account [21]. At high temperatures, the right-hand side of Eq. (4) 
dominates the Hubble expansion term and the WIMP number density is forced to follow 
its equilibrium solution, n ^ ncq, which for T < is well approximated by a Maxwell- 
Boltzmann distribution. As T drops further, however, the Boltzmann suppression of n 
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(and neq) means that the WIMP annihilation (and creation) rate will eventually drop 
below the Hubble expansion rate; at chemical decoupling, T = T^, the co moving number 
density a^n starts to deviate considerably from its equilibrium value and finally freezes 
out and stays constant; for WIMPs, this happens around T^d ~ 

Even at temperatures below Ted, however, the WIMPs are kept in local thermal 
equilibrium by scattering processes with SM particles (for which the full collision term is 
provided in Appendix A). As pointed out in [22], it is sufficient for the discussion of the 
kinetic decoupling process to consider the second moment of the Boltzmann equation 
(instead of the first moment like for chemical decoupling): multiplying Eq.(3) by p'^/E, 
integrating it over p and keeping only the leading order terms in p^/m^ results in 

id, + 5H)T^ = 2m^ciT)iT-T^) , (5) 

where c(T) is given in Eq. (A. 8) and is defined by 

^^pV(p) = 3 m^T^n^ . (6) 

Note that this "temperature parameter" for the WIMPs does not require any 
assumptions about the form of /(p), but simply provides a convenient means of 
characterizing the deviation from thermal equilibrium (for which T^ = T holds). 
During radiation domination, the expansion rate is given by 
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where Qes is the effective number of relativistic degrees of freedom. Introducing further 

x = m^/T, (8) 

y ^ m^g-J"TjT' , (9) 

one can now bring Eq. (5) into the convenient form 



with g^l"^ = glg^/ll + \^T). From this equation, one can directly read off the 
asymptotic behaviour of the WIMP temperature: At large T, the term in front of the 
right-hand side is much larger than unity and one has = T; when T becomes small, 
the WIMPs finally decouple completely from the thermal bath and y stays constant, i.e. 

oc T'^glf^ oc a~^. The kinetic decoupling temperature is thus conveniently defined by 
equating these two regimes, as if the decoupling process were to occur instantaneously: 

Xkd = 7^ = ^cfr(Tkd) • (11) 

This procedure is illustrated in the left panel of Fig. 1 which shows, for a typical WIMP 
candidate, the phaseplot and solution of Eq. (10). As one can see, kinetic decoupling is 
indeed a process that happens on a rather short time scale; it becomes also apparent that 
thermal equilibrium is maintained extremely efficiently for T < T^d, with = T being 
an attractor solution that would be restored almost immediately after any departure 
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Figure 1. The left panel shows the phaseplot and solution for the WIMP temperature 
evolution, for ^ 100 GeV and \M\ ^ (^y (to^/cj)^, expressed in the dimensionless 
variables introduced in Eqs. (8,9). At T < Tkd, any departure from thermal 
equilibrium (T^ = T) is restored almost immediately (except for a short period around 
the QCD phase transition); for T > Tkd, the WIMPs decouple from the thermal bath 
and cool down with the Hubble expansion as cx a~^. 

In the right panel, the effective number of rclativistic degrees of freedom is plotted 
as a function of the temperature, implementing the results of [25] for the evolution of 
this quantity during the QCD phase transition; for reference, the decoupling of muons 
and electrons is also indicated. 



from this behaviour (except for a short period during the QCD transition, see below, 
when the rapidly changing effective number of degrees of freedom does not allow this). 

In principle, the scattering with all types of SM particles contributes to c(T), see 
Eq. (A. 8). This picture is a bit complicated by the fact that kinetic decoupling in some 
cases can take place close to, or even above the QCD phase transition, the details of 
which are not yet fully understood. Lattice calculations, however, start to converge at 
a value for the critical temperature of ~ 170 MeV for the most interesting case of 
two light (up and down) and one more massive (strange) quark flavour [23] and indicate 
that the plasma can be described by free quarks and gluons only for T > 4Tc [24]. For 
the effective number of degrees of freedom during the transition, we adopt the results 
of [25] as displayed in the right panel of Fig. 1. As scattering partners are concerned, 
we conservatively restrict ourselves to leptons and, for T > 4Tc, to the three lightest 
quarks. 

The resulting range in T^d for neutralino dark matter, obtained after having 
performed the extensive scan described in Section 2, is shown in Fig. 2 as a function of 
the mass and gaugino fraction Zg = lA^^nl^ + |A'^i2|^ (in our case dominated by the 
Bino fraction). The gray band indicates the QCD phase transition; values for T^^ inside 
or above this band should be interpreted as upper bounds on the decoupling temperature 
since the scattering with some of the hadronic degrees of freedom was not taken into 
account. On the other hand, as the coupling of WIMPs to hadrons is usually smaller 
than to leptons, the difference between this upper bound and the actual value of T^d is 
not expected to be very big; note also that the scattering with bound QCD states like. 
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Figure 2. The range of decoupling temperatures for ncutralino DM. For models that 
fall inside or above the gray band marking the QCD phase transition, the actual value 
of Tkd will be slightly smaller than indicated. See text for further details. 



e.g., pions is suppressed due to their rather large masses and thus small abundance (the 
evolution of density fluctuations, on the other hand, may very well be influenced by the 
details of the QCD phase transition, see the next Section). In addition to the result 
of the scan, the Figure also indicates the decoupling temperature for four mSUGRA 
benchmark models that were introduced in [26, 16] and present typical examples for 
neutralinos in the bulk (/'), coannihilation (J*), funne\{K') and focus point (F*) region; 
the quite different annihilation spectra in gamma rays for these models, and the resulting 
prospects for indirect detection, were recently studied in some detail in [16, 27]. 

Assuming a constant equation of state and relativistic scattering partners, Eq. (5) 
has actually an analytic solution [22] that can be used for a quick estimate of Tkd- 
This estimate proves to be rather good (within 10% of the full result shown here) for 
masses below a few hundred GeV; above that, however, the exact mass dependence of 
the number density of r leptons, in particular, can be crucial, leading in some cases to 
differences of more than a factor of 5 between the two results. A further observation 
is that MSSM and mSUGRA models occupy roughly the same regions in the Tkd ^ 
plane; the largest values of Tkd for Higgsinos (at high masses) and Binos (at intermediate 
masses), however, corresponds nearly exclusively to MSSM models, while almost all 
Binos with a decoupling temperature below the band occupied by mixed neutralino are 
mSUGRA models. 

For illustrative reasons, flnally, the right panel of Fig. 2 compares the kinetic with 
the chemical decoupling temperature: as anticipated, kinetic decoupling takes place 
much later than chemical decoupling, at temperatures a factor of 10 — 1000 lower. The 
possible range in Tkd is, furthermore, considerably larger than the one in Ted - which 
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of course simply reflects the fact that the DM relic density is constrained extremely 
well while there are so far no observations that would put stringent bounds on (see 
also Section 5.1). Note also that Tkd is fairly uncorrelated with the chemical decoupling 
temperature, as well as with the neutralino annihilation cross section. 

To conclude this Section, let us recall that the formalism presented here keeps 
the leading order terms in p^/m^, thus allowing the determination of the decoupling 
scale to an accuracy of O (x^^); while this is usually more than sufficient, it would be 
straightforward to include also higher orders in the discussion of the Boltzmann equation 
and thereby arrive at an even more accurate description of the decoupling process. 



4. From decoupling to the first protohalos 

Before kinetic decoupling, small-scale perturbations in the CDM fluid are damped due 
to the tight coupling of the WIMPs to the heat bath. For temperatures T < Tkd, flrst a 
remaining viscous coupling between these two fluids and then the free-streaming of the 
WIMPs generate an exponential cutoff in the power spectrum [13], with a characteristic 
comoving damping scale of 



\ 1/2 / 

TmJ ln(4aeq/akd) ao ""^ 



He^. (12) 



The WIMP mass contained in a sphere of radius n/kfs is thus given by 
4n /ttV J ^ + (4'^kd/50 Mev) /19.1 



Mfe ^ —p.y — = 2.9 X 10"^ ^ —-7^ ^ — Mq . (13) 

Later, it was shown that acoustic oscillations also have to be taken into account [14, 15]. 
Leading to a similar exponential cutoff, the characteristic damping mass in this case is 
given by the CDM mass inside the horizon at the time of kinetic decoupling: 

/ 1 /. \ -3 

Mao - 3 ^3 



3.4xlO^M^^^l Mq. (14) 



T=Tkd \ 50 MeV 

Note that the above expressions contain the full dependence on Qcs (evaluated at Tkd), 
which is often not quoted despite the large range of possible decoupling temperatures; 
to arrive at the numerical values, the most recent estimates for fl^h"^ and Qmh"^ [1] were 
used. 

The two cutoff scales are plotted against each other in the left panel of Fig. 3. As 
noted before in [14, 15], the effect of acoustic oscillations is, indeed, more important for 
the case of a Bino-like neutralino with Xkd ^ 10^ (these models occupy the top part of the 
Figure, just to the right of the band containing mixed neutralinos) . In general, however, 
this statement turns out to depend on the WIMP nature and there are cases where the 
free-streaming cutoff is almost 2 orders of magnitude above the one induced by acoustic 
oscillations. Since the two damping mechanisms are physically independent, the actual 
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Figure 3. The left panel shows the exponential cutoff scales associated to the main 
damping mechanisms of the matter power spectrum after kinetic decoupling, viz. free 
streaming and the effect of acoustic oscillations, respectively; for models above (below) 
the dashed line, the former (latter) mechanism thus provides a stronger suppression of 
the power spectrum. In the right panel, the cutoff mass resulting from the dominating 
of these two independent effects is plotted against the neutralino mass, indicating the 
typical size of the smallest protohalos to be formed. 



cutoff Mcut in tlie power spectrum is tlius, ratlier, given by Mcut = max [Mfs, Mao]; tlie 
possible range of Mcut is displayed in tlie right panel of Fig. 3 as a function of m^. 

For very small values of Mcut, corresponding to large T^d, one might wonder whether 
the QCD transition could leave an imprint on the power spectrum. In fact, if it is first 
order, the sound speed vanishes during the transition and density perturbations fall 
freely, potentially leading to the production of DM clumps with masses of 10~^° to 
10~^°Mq [8]. However, the corresponding enhancement factor in the CDM density 
fluctuations is only between 2 (from a lattice fit) and 20 (using the bag model) at scales 
of ~ 10~^^Mq and significantly smaller at larger scales; this has to be compared to the 
exponential suppression of power below Mcut due to the damping mechanisms discussed 
here. For the smallest cutoff scales shown in Fig. 3, Mcut ^ IO^^'^Mq, the actual cutoff 
mass might thus be slightly, but certainly not very much, smaller than indicated - 
depending on the details of the QCD phase transition. 

Following the paradigm of hierarchical structure formation, the smallest scales, 
and thus the scales closest to the cutoff, typically enter the non-linear regime first. The 
smallest gravitationally bound objects to be formed in the universe are in that case also 
the first; protohalos with a mass of around Mcut- This behaviour has been confirmed 
numerically, where these protohalos could be followed until a redshift of ^ ~ 26 [28]. 
The range of expected minimal protohalo masses displayed in Fig. 3 is only slightly 
smaller than what was found earlier [29] using an order-of-magnitude estimate for Tkd 
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(based on [13]) instead of the exact value as defined by the solution of Eq. (10). For 
a given model, however, it turns out that the difference in the infered cutoff mass still 
is typically about a factor of 10, rather independent of m^; adding to this the effect 
of identifying Mao (like in [14, 29]) instead of M^nt = max [Mfg, Mao] with the smallest 
protohalo mass, this difference can in some cases increase to a factor of almost 1000. 

5. Discussion 

5.1. Observational prospects 

Observational prospects depend crucially on whether the first protohalos survive until 
today or whether they are disrupted due to tidal interactions in merger processes or 
encounters with stars - an issue that is still under debate. Several studies show that 
even if the protohalos lose some of their material on their way, most of the mass resides in 
a dense and compact core that remains intact [28, 30]. Other studies are less optimistic 
[31]. However, one should keep in mind that the very first objects probably form out of 
rare fluctuations with Mrf > Mcut that enter the non-linear regime before fluctuations 
of size Mcut; these DM clumps could develop denser cores before serious encounters, 
enhancing thus their survival probability. The importance of rare fluctuations for the 
formation of microhalos is also supported by numerical simulations [32]. 

A clumpy halo has important consequences for the indirect detection of DM, where 
one tries to discriminate exotic, DM-induced contributions to cosmic rays of various 
kinds from the standard astrophysical background. DM annihilating in dense clumps 
might, in particular, be visible in the form of gamma-ray point sources. In general, 
however, one can only expect to see some of the larger clumps [33] with gamma-ray 
telescopes like Fermi/GLAST [34], while the very small clumps with masses close to Mcut 
are unlikely to be resolved [35]. They do, however, contribute to the diffuse flux; indeed, 
it has been argued that models where Fermi is expected to see some of the larger clumps 
generally feature so many small clumps that their contribution to the diffuse flux would 
be in conflict with already existing observations [36]. In any case does the so-called boost 
factor, i.e. the enhancement in the smooth component in a clumpy halo with respect to 
the flux expected for a homogenous DM distribution, depend on the small-scale cutoff 
of the clump-distribution - which becomes even more important when one tries to fully 
take into account the self-similar nesting of halos (subhalos within subhalos). The actual 
functional dependence on Mcut, however, strongly depends on the density profile of the 
smallest subhalos today as well as on their radial distribution in the Milky Way. Since 
the resolution of even the so far most ambitious numerical A^-body simulations is still 
way above the expected cutoff scale, one has to extrapolate these results by many orders 
of magnitude, the implications of which are, understandably, still under debate. The 
most likely dependence of the boost factor on Mcut is logarithmic, i.e. each decade in 
halo masses contributes about the same to the total boost. The boost factor for charged 
cosmic rays, conceptually a bit different from the one for gamma rays since charged 
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particle propagation through the diffusive halo has to be taken into account, may be 
much less than for gamma rays [37]. 

It has been proposed that the proper motion of the smallest subhalos could be 
detectable [38], but prospects for this idea do not look very promising even for the 
most favourable case of a rather large Mcnt ~ 10~'^Mq; again, it is the EGRET diffuse 
background that places severe restrictions on this idea [39]. While direct detection 
rates for DM would of course be greatly enhanced if a DM clump happens to pass 
through the earth, prospects are on average worse than for a smooth DM distribution 
[40]. Microhalos close to the cutoff scale, furthermore, feature a virial radius much 
larger than their Einstein radius, so gravitational lensing is commonly not considered 
an adequate method of assessing these smallest structures, either; recently, however, it 
has been suggested that a detailed study of multiple images of time-variable compact 
sources in strong lensing systems actually may have the potential to constrain the power 
spectrum in the lensing system even at scales close to typical values for the cutoff [41]. 
The maybe most promising way to directly probe the small-scale distribution of DM 
might be to use anisotropy probes such as the angular correlation in the extra-galactic 
gamma-ray background [42] or the gamma-ray flux (one-point) probability function [43]. 
While galactic substructures do provide the dominant contribution to the first signal, 
they do so mostly for rather large clumps with masses M > IO^Mq [44], which would 
imply that this method cannot be used to asses the small halos we are interested in 
here. The very recently proposed second method, however, claims to probe the subhalo 
distribution, at least in principle, down to the smallest scales. 

To conclude this Section, let me mention that the role of microhalos in structure 
formation, and the implication of their presence for dark matter experiments, is also the 
subject of another review article in the present NJP focus issue on dark matter [45]. 

5.2. N on-neutralino dark matter 

All expressions in Sections 3, 4 and Appendix A are equally applicable to any WIMP 
candidate and the resulting decoupling temperatures and cutoff masses are expected to 
always fall into the range presented in the preceeding sections. In fact, in many cases 
the range will be much smaller simply because there are less free parameters to play 
with than in the case of supersymmetry. As an example. Fig. 4 shows the situation 
for Kaluza-Klein DM in universal extra dimensions [6], where the mass splittings 
between the Kaluza-Klein particles arise due to a combination of electroweak and 
radiative contributions and are generally expected to be rather small. In the minimal 
setup [46], the whole spectrum of states can be specified by only two parameters; the 
compactification scale R and the cutoff scale A beyond which the four- dimensional 
effective theory ceases to be valid. For Kaluza-Klein dark matter, the cutoff in the 
matter power spectrum turns out to always be dominated by acoustic oscillations and, as 
can be seen in the Figure, falls into a very narrow range around Mcut = Mao ~ 1O~^M0. 
Another interesting class of DM candidates are particles that, unlike WIMPs, have 
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Figure 4. The kinetic decoupling temperature and the corresponding cutoff scale 
for the lightest Kaluza-Klein particle in universal extra dimensions. For high Higgs 
masses, rrih > 150 GeV, the grey region corresponding to the relic density constraint 
given in Eq.(2) shifts upwards, allowing an LKP mass up to around 1 TeV [47]. 



never been in thermal equilibrium and for which the formalism presented here therefore 
does not apply (though it would be very interesting to extend it such as to cover even 
these cases). The maybe best studied example is the axion [48]. Initially misaligned, it 
starts to oscillate coherently around T ~ 1 GeV and from then on evolves as pa oc 
just like ordinary CDM. Inhomogeneities on scales larger than the Hubble horizon at 
the temperature of realignment later evolve into axion mini-clusters with a typical mass 
around 10~^^Mq [49]. One should, however, keep in mind that for particles that decouple 
before the QCD transition - like axions, primordial black holes with Mpbh ^ Mq 
or ultra-cold WIMPs [50] - small-scale fluctuations may be strongly amplified if the 
transition is first order, producing DM clumps with masses of 10"^° to 10^^^ Mq [8]. 

Finally, if DM consists of superWIMPs that result from the late decay of thermally 
produced WIMPs, the actual cutoff in the power-spectrum is not the one from the 
WIMP decoupling but the one that is imposed from the kinematics of the decay (through 
the mass difference between decaying particle and DM particle). In fact, such models 
have been proposed to address a certain tension that is sometimes claimed at "small" 
scales (in this case Mpc instead of the pc scales that correspond to Mcut ~ 10~^Mq) 
between observations and numerical iV-body simulations [51]. However, this idea works 
only partially [52]; what is more, the evidence for small-scale "problems" of standard 
ACDM cosmology may soon well disappear completely, with more detailed observations 
and A^-body simulations starting to converge [53]. Nevertheless, late-decaying DM is 
an interesting possibility that does not have to be related to this particular idea; in 
contrast to the typical Mcut for WIMPs, a large cutoff in the power spectrum might 
even be possible to probe by future micro-lensing missions. 
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6. Conclusions 

The kinetic decoupling process of WlMPs from the thermal bath can be followed in great 
detail by solving the full Boltzmann equation in this regime. Extending the formalism 
presented in [22], by allowing for non-relativistic scattering partners and taking into 
account the full time-dependence of the effective number of degrees of freedom, a highly 
precise determination of the decoupling temperature becomes possible that in turn can 
be translated into a small-scale cutoff in the spectrum of matter density fluctuations. 

An extensive scan over the parameter space for supersymmetric neutralino DM 
reveals a slightly smaller range of cutoff masses, 10~^^Mq to a few times 10~'^Mq, but 
basically confirms the only existing corresponding scan so far [29], which is based on 
an order-of- magnitude estimate for the decoupling temperature (given in [13]). The 
resulting difference in Mcnt for individual WIMP models, however, can be sizable; 
typically of the order of 10, models with a difference of almost 10'^ were found. Another 
important result of the scan presented here is that whether free streaming or acoustic 
oscillations are more effective in the suppression of power on small scales depends on 
the DM particle nature (in slight disagreement to the claim of [15] who presented a very 
detailed study of the evolution of density contrasts through and after kinetic decoupling, 
albeit based only on one particular DM candidate). 

The range of decoupling temperatures and cutoff masses presented here is indicative 
for the whole class of WIMP DM candidates, though many models - such as Kaluza- 
Klein DM - will exhibit a much smaller range. For non-WIMP candidates, the mass 
of the smallest clumps can differ significantly from the range derived here; it would be 
interesting to develop tools that allow an as precise determination of the cutoff scale 
for these CcLSGS cLS for the case of WIMPs. As for detectional prospects of the smallest 
DM clumps, many interesting ideas have been put forward. Though challenging, it is an 
exciting possibility that one may be able to measure the DM distribution on such scales 
in the future. In order to really address the connection to the microphysics of the DM 
particles, however, one still needs a better understanding of how the first protohalos 
evolve and, given their initial distribution, what they are expected to look like today. 

The routines for calculating the kinetic decoupling temperature and the associated 
cutoff scale have been implemented in DarkSUSY [17] and will be available with the next 
release [20]. 

Appendix A. The collision term in the Boltzmann equation 

In this Appendix, we derive the collision term for scattering processes between non- 
relativistic WIMPs X considerably less massive (usually SM) particles that are in 
thermal equilibrium with the plasma in the early universe. For the former, we use 
p'^ = {E,p) to denote ingoing momenta, while for the latter we use k'^ = (co',k); the 
corresponding outgoing quantities are marked with a tilde. The treatment here follows 
closely that of Ref . [22] , but slightly extends it in taking into account the possibility of 
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non-relativistic scattering partners. 

The collision term, i.e. the right-hand side of the Boltzmann equation, reads 



where 



{2Tif2uj J {2Tif2Co J {2'Kf2E 

J ^gsM[[lT g^iuj)) g^{^)f{p) -{IT g^{^)) ^7^(^)/(p)] , (A.2) 
gH-^)= {e-'^±l)-\ (A.3) 



Here, \M.\ is the scattering amplitude squared, summed over final and averaged 
over initial spin states, and g^m is the number of internal degrees of freedom of the 
scattering partner (the upper signs apply for fermions, the lower for bosons); as for the 
following expressions, a summation over all SM particles in thermal equilibrium is always 
understood. No assumptions about the x distribution function /(p) are necessary; as 
long as the WIMPs are much less abundant than their scattering partners, however, 
Pauli suppression factors for / can safely be neglected - as has been done in Eq.(A.2). 

After chemical freezeout, one typically has ~ T ^ m-y^. For kinematical reasons, 
the average momentum transferred during the scattering events is thus small, so (A.l) 
can be expanded as 

CXD 

C(E)=^C^ (A.4) 
d^k f d^k f d^p 



(27r)32cj7 (27r)32cD7 {2tt)32E 



■2 



X {27r)^6{E T Co - E - uj)\M\ J 



iz)^(p)5(3)(p-p) 



(A.5) 



where 



Dq(p) = q.Vp = (k-k)-Vp, (A.6) 

and the derivatives of the delta function are as usual defined through an integration by 
parts. Following similar steps as detailed in Ref. [22], but paying special attention to 
the difference between u and k = |k|, one can now proceed to calculate the expansion 
coefficients . To lowest non-vanishing order in p'^/E'^ and u/m^, one then arrives at 
the following expression for the collision integral: 

C = C' + C^ = c{T)ml [m^TAp + p ■ Vp + s] /(p) , (A.7) 



where 



^^^^ ^ 6(27r)3m4T 



^ ^ I dkk^'u-^g^ {It g^)\M\]^, . (A.8) 

For clarity, the sum over all SM scattering partners i has here been made explicit. 

As an interesting aside, the above integral can be solved analytically in the 
relativistic limit, mi 0, if the amplitude squared scales hke a power of u [22]. For 



Particle DM and small-scale structure 



15 



WIMP scattering below s-channel resonances, ^ ^ cjrcs, e.g., the scattering 

2 2 

amplitude is always given by \ Ai\ = \.M\q (uj/m^)'^, see Appendix B. In this particular 
case, the integral in Eq. (A. 8) becomes: 

J dk k'uj-' (1 T 9^) Wf "^^<^"-= 720 C(7) Mo ^x' (^-9) 

for bosonic scattering partners; for fermions, the above expression has to be multiplied 
by 31/32. 

Appendix B. Neutralino-fermion scattering amplitude 

The scattering between neutralinos and fermions is mediated by sfermion exchange in 
the s- and w-channel, and Z- and Higgs-boson exchange in the t-channel. Denoting with 
p (p) the ingoing (outgoing) momentum of the lightest neutralino Xi? and with k (k) the 
ingoing (outgoing) momentum of the fermionic scattering partner i, the total scattering 
amplitude is given by 

M = Ms + Mu + Mt, (B.l) 

with 

Ms = ^J2 ulC-^ {^f;P^ + 4;P«} % {4^P, + 4^P^} Cul (B.2) 

i=l,2 
i=l,2 

Mt = iAt{Z) Upgznl57f^UpU-^Y {gzee^L + gzeePn} 

-i At{Hs) Up {3? [gn^u] - 275 ^> bLn] } ^P^gH^u^k 

m=l,2 

+iAt{H3) Up {3ft [gH.,11] - [gHsu] } Up U-^gH3Ul5Uk , (B.4) 

and 

A(s,t,u){A) = - [{s, t, u) -m\ + im^TA] ~^ (B.5) 

for any particle or sparticle A with mass rriA and width Ta- In the above expressions. 
Hi and H2 denote the scalar Higgs-bosons and its pseudoscalar version, C the 
charge conjugation matrix and a superscript T the transpose of a matrix. For the 
coupling constants, the same conventions as in DarkSUSY [17] are adopted (for an explicit 
representation of these couplings, see e.g. [54]; for a nice summary of Feynman rules for 
Majorana fermions, see [55]). In the case of neutrino scattering, there is only one 
sneutrino to be exchanged; one thus has to drop the corresponding sums over sfermion 
mass eigenstates (i, j) and to replace iij —^uin the above and the following expressions 
(in this case, one may of course also safely neglect the exchange of Higgs bosons). 

The typical CMS fermion energy is much smaller than the neutralino mass, 
i.e. u <^ m^. As explained in Appendix A, we are therefore only interested in small 
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momentum transfer, 



t^O, (B.6) 
s — > + 2m^ijj + ml . (B.7) 
In this limit, averaging over initial and summing over final state spins gives: 

= + + + 23fJ [MMt - MsMl - MtMl] , (B.8) 

where 

(B.9) 

9 9 

(B.IO) 
(B.ll) 



|A^„| = \Ms\ {to -u) 



i,m=l,2 



,2 



2'^ 



m} 

1 



X 



9^ 



53 



9Lp^ 9 hp 



9iiti9i^a 



at 



2 2cj^ 
+ — sft 



ill 



1=1,2 



'^9z£e 



9hp^ 9 hp 



9p,a^ - 9i^a'^i 



^9z 



Zll 



(B.12) 

2 



9zpp 



12 2 

[uj — mg 



i=l,2 m=l,2 

-25R L^:^? 







2 


















9i^p\ 


+ 


^iip\ 





mpijj 



iiti^tiPi 



mo 



(B.13) 

MtMl = -MtMliiu -iu) (B.14) 

The above expressions extend the formulae given in [9] for neutralino-neutrino scattering 
and correctly reproduce the expressions obtained in [56] for neutralino annihilation (after 
the replacements t u,u ^ s, s ^ u). 

For relativistic fermions, the above expression further simplifies considerably. Below 
the resonance(s) 



^res,i — 



2 2 

m| — mr. 



2m^ 



the scattering amplitude becomes 



l-MI 



IjJ 



m^ 



m 



X 

mz ) 



Y 



^7lii 



9^11 + fl'fa 



(B.15) 



(B.16) 
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m?- /rn^ 



Above the resonance, c^res ^ ^ itl^, the s- and -u-channel amphtudes are roughly 
constant, while the t-channel amplitude continues to scale like Ait oc u. In the limit 



of a pure Bino (i.e. 
sfermions, Eq.(B.16) 



V2gYYs., 



V2gYY, 



diQzil 



Decomes the expression earlier found in [10]: 



\M\ 



■2 (Bino) 



UJ 



m- 



0), and degenerate 



(B.17) 



Let me stress that the expressions (B.16) and (B.17) are just given here for convenience; 
for all numerical calculations, the full expression (B.8) was used. 
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